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Abstract 

We assess the potential of detecting a charged Higgs boson of the MSSM at the LHC 
via its decays into a chargino and a neutralino. We focus our attention on the region of 
parameter space with m H ± > mt and 3 < tan/3 % 10, where identification of the via 
other decay modes has proven to be ineffective. Searching for means to plug this hole, 
we simulate the decays — > xfxi an d — ► XiX2>XiX3 — the former can yield 
a single hard lepton (from the chargino decay) while the latter can yield three leptons 
(from the chargino and neutralino decays). Coupled with the dominant top quark + 
charged Higgs boson production mode, the resulting signature is one or three hard, 
isolated leptons, substantial missing transverse momentum and a reconstructed (via a 
3-jet invariant mass) top quark. The single lepton channel is swamped by background 
processes; however, with suitable cuts, a trilepton signal emerges. While such a signal 
suffers from a low number of surviving events (after cuts) and is dependent on several 
MSSM input parameters (notably M 2 , u, and slepton masses), it does fill at least some 
of the void left by previous investigations. 
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A pair of spin-0 charged Higgs bosons, H ± , arises in any Two Higgs Doublet Model 
(2HDM) alongside a trio of neutral Higgs bosons — the CP = +1 'light' h and 'heavy' H 
(with rrih < uih) and the CP = — 1 'pseudoscalar' A. The charged Higgs bosons have been 
at the focal point of extensive studies since they have no Standard Model (SM) counterpart, 
and thus could provide irrefutably clear evidence of an extended Higgs sector and new 
physics beyond the SM. On the other hand, it may be difficult to either distinguish one 
type of neutral 2HDM Higgs boson from the SM Higgs boson or observe more than one 
of the neutral species [|1], 0, 0, [|. Embedding the 2HDM inside the attractive theoretical 
framework of supersymmetry (SUSY) yields (together with additional assumptions about 
minimal field content and minimal number of new couplings) the Minimal Supersymmetric 
Standard Model (MSSM). In the MSSM, at tree-level, the masses of all Higgs bosons, along 
with their couplings to the SM fermions and gauge bosons, can be parametrized in terms of 
only two unknown input parameters, generally taken as the mass of one of the Higgs bosons 
(typically either rriA or m H ±) and the ratio of the vacuum expectation values (VEV's) of the 
up-type and down-type Higgs doublets (denoted by tan/5) 0. 

As is well-known M, these tree- level relations can receive substantial radiative correc- 
tions, most importantly enabling rrih > Mz (making the upper limit on rrih in the MSSM 
~135GeV 0). However, the tree-level relation between the masses of the charged Higgs 
bosons and A, m 2 H ± = m\ + M^ v± , is almost invariably quite insensitive to such correc- 
tions f8j. Properly taking into account the corrections to the light Higgs boson mass, m H ± 
may still be related to rrih, and an indirect lower bound of ~140 GeV || for tan/5 ~ 3-4 can 
be placed on the former due to the thus-far null search for a Higgs boson at LEP2. This 
bound grows rapidly stronger as tan/3 is decreased while tapering very gradually as tan/? is 
increased (staying in the 110-125 GeV interval for tan/3^5). There are also other processes 
where charged Higgs bosons (or A , to whose mass that of the H ± is closely tied) enter as 
virtual particles at the one-loop level. These include neutral meson mixing (K°K°, D°D°, 
or B°B°) |I| Z° -> bb (R b ) P, ||], and b -> 57 decays |I| 0, U, The b -> 57 



decays are generally thought to be the most constraining ||T2[ . Here restrictions on m H ± are 
linked to a number of MSSM variables, notably including the masses of the lighter chargino 
and the stops. We reserve further comment on this somewhat complex issue until after we 
present our results. 

More direct limits on charged Higgs bosons come from hadron collider searchesQfor lepton 
non-universality (excess taus) in top quark decays resulting from t — > bH + followed by 
H + — > t + u t |TBj and the charge-conjugate reactions^ (excessive numbers of charmed final 



states in top decays resulting from H + — > sc may also be of use, as well as H + — > W + bb |L7| 



for tan/5 ~ 1). At the soon to commence Run II of the upgraded Fermilab Tevatron, such 
channels will allow experimenters to scan the MSSM parameter space for large and small 
values of tan/3 roughly up to the kinematical limit of the t — > bH + decay, m t — [|18[]. The 



reason for the tan/3 dependence stems from the couplings between a charged Higgs boson 



1 There are also direct searches for charged Higgs boson pair production at LEP2; however, bounds 
obtained in this way are relatively low, m H ± > 77.5-78.6 GeV [Q. At a future 500 GeV e + e~ linear collider, 
this could increase to a potentially competitive ~210GeV ]l5| ]. 

2 Hereafter, inclusion of charged-conjugate processes may be assumed unless explicitly excluded. 



1 



and top and bottom quarks, given by[] 



g 2 

~ , f 2 H + (m t cot(3tb L + m b t&n(3tb R ) . (1) 



The square of this, which has a minimum at tan/3 ~ 6 -7, is proportional to the strength of 
either the gb — > iH + cross-section or the t —>■ bH + decay width. In the intermediate tan/? 
region around this minimum, a Tevatron search for charged Higgs boson pair production, 
qq — > H + H~ , which mainly proceeds utilizing only gauge couplings, could be feasible fl9j if 



the charged Higgs bosons are light enough — certainly m#± <m t . As the mass of the charged 
Higgs boson grows larger than m t , simple phase space suppression will severely handicap pair 
production. 

Thus the likely legacy bequeathed to the Large Hadron Collider (LHC) will be the pursuit 
of a heavy charged Higgs boson (with m H ± >m t ). At the LHC, the dominant production 
mechanism for heavy H ± scalars is via the 2^2 reaction gb — > tH~ [^] and the 2 — > 3 
reaction gg — > tbH~ [^] . Alternative production modesf] are charged Higgs pair production, 
gg,qq — > H + H~ f2"3| , and associated production, gg,qq — > W ± H T f2"5fl . The former suffers 
from a lack of phase space for m H ± > m t and low quark parton luminosities (qq) in the LHC 
protons or heavy propagator loop suppression (gg). The latter suffers from a huge irreducible 
background induced by either ti and/or W + W~ production and decay, depending upon 
whether the charged Higgs boson decays via bi, W~h or t~v t states [^]. For the preferred 
production mechanism, henceforth to be collectively referred to as the top-H^ production 
mode, the tan/3 dependence follows from the square of Eqn. (HD, again making high and low 
tan/3 values more accessible. The connection between the 2 — > 2 and the 2^3 reactions 



has been discussed numerous times before 27, 2S|: the former is obtained from the latter 



if one of the gluons splits into a bb pair, with one bottom quark (i.e., b- or 6-quark) then 
interacting with the remaining gluon while the other is assumed to act as a spectator. The 
appropriate procedure |28|, |2P| for estimating the inclusive Higgs production cross-section is 
to combine the 2^3 subprocess with the 2 — > 2 subprocess through the subtraction of a 
common logarithmic term ~ a s \og(Q 2 /m%). However, utilizing the 2 — > 2 simulation for the 
kinematical event selection (and thereby tacitly assuming that the final state bottom quark 
manifest in the 2^3 formulation is soft and thus untaggable) could lead to erroneous 
event-shape parameter distributions (transverse momentum, opening angles between jets, 
etc.) due to the possible presence of the extra, neglected 6-jet. Therefore, here we simulate 
both subprocesses: 2 — > 2 simulations are employed solely to normalize the cross-sections, 
making use of this subtraction procedure, while 2^3 simulations are used to implement 
our selection and acceptance cuts with the 6-jet resulting from the bottom quark produced 
in the 2 — > 3 reaction subject to the same acceptance, resolution and isolation criteria as 
any other jet in that event. 

Several decay modes for heavy states (for branching ratio studies, see |30[), have been 



analyzed assuming the above top-if production mechanism, including: H — ► bt [27, 31 



3 Analogous formulae hold for the other two SM fermion generations. 

4 The bq — > bH^q' mode of can only be relevant for very high values of tan/3. 
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generally expected to be the dominant decay mode; H — > sc (in ATLAS study of ||31|| ); 
H — > VK~/i f32|l ; and — > 7"~z/ T f33| . All these decay modes (that will hereafter be termed 
"SM" decays) were simulated (including parton shower, hadronization and detector effects) 
in either the ATLAS simulations of pT] , p2| . |33| and/or in f2"0fl , with the latter concluding that 
H^ 1 scalars with masses up to ~400 GeV can be discovered by the LHC, but only if tan/3 < 3 
(which is in the neighborhood of the indirect limit from LEP2) or tan/3 > 25. The ATLAS 
studies roughly concur, adding that the H — > t~v t channel can push the high tan/3 reach 
down below 20 for m#± <400 GeV, with a minimum at tan (3 ~ 10 when m#± is close to m t . 

The purpose of this paper is to assess the prospects for utilizing the thus-far neglected 
SUSY decay channels of the MSSM charged Higgs bosons to probe regions of the parameter 
space inaccessible at the LHC via the SM decay modes. The charged Higgs boson can 
in fact decay predominantly via these SUSY modes, as illustrated in Fig. 1 — which also 
serves to highlight the potential significance of these SUSY channels in the large m H ± and 
intermediate tan/3 regions. In particular, we explore H ± decays into a chargino (xt) an d a 
neutralino (xj)'i i- e -i H ± — > xtx°ji i = 1 or 2 and j = 1, 2, 3, or 4. The lightest neutralino, 
Xi, is assumed to be the stable lightest SUSY particle (LSP). The decay width is given by 



35fl: 



- xt%) = m X ~ 3 J (2) 

where F L and F R are as follows: 

F L = cosPiNjtVa + ^iNjz + N^ttme^Vx], 

F R = sin^N^Un- ^ 2 (N j2 + N jX tan 6 w )U i2 ]. (3) 



(For the U, V and iV matrices, we have followed the notation of 1331.) Here, m~±(m~o) are 



the masses corresponding to the xt(Xj) states and e is the sign convention for the neutralino 
mass eigenstates. Dependence on the additional MSSM input parameters M 1) M 2 , and // 
enters from the gaugino/Higgsino mixing matrices via the U, V and N. M 1 and M 2 are 
the U(1) Y and SU(2) L gaugino masses, respectively, and /x is the Higgsino mass parameter. 
Grand Unified Theories (GUT's) predict gaugino unification and M 1 = | tan 2 9 w M 2 , as will 
be assumed in all numerical calculations. 

Branching ratios (BR's) for the chargino-neutralino decay modes of the charged Higgs 
bosons are shown (along with the important SM decay BR's) versus tan/3 in Fig. 2, choosing 
M 2 = 200 GeV and /j = —120 GeV as in Fig. 1. While this point is favorable for chargino- 
neutralino decays, it did not result from a exhaustive search for the optimal choice. Three 
charged Higgs boson masses are examined {m H ± = 200, 300, and 400 GeV). For m H ± = 
200 GeV, the only chargino-neutralino decay channel open is XiXv whereas, for m#± = 
300 GeV, the X1X2 an d X1X3 channels are also accessible. In fact, in this latter case the 
BR for H ± — > X1X2 1S larger than that for H ± — > X1X1 f° r tan/3 £2. By the time m#± 
reaches 400 GeV, many chargino-neutralino decay channels have opened up and the situation 
becomes fairly complicated. Decays to the heaviest charginos and neutralinos may well 
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generate cascade decays rather than (predominantly) decaying directly to the LSP. This will 
introduce additional MSSM parameter space dependence as well as complicate the analysis. 
Further, as we shall see, for if ± masses much beyond this point, the lower top- H ± production 
rate robs us (after the necessary cuts) of any signal events in the multilepton channels we 
will be investigating. Therefore, there is considerable justification for concentrating upon 
the H ± — > xt Xi) X1X21 XiXjs channels in this exploratory study. 

Note from the m^± = 400 GeV plot in Fig. 2 that the sum of the various chargino- 
neutralino modes (which is represented by the "all SUSY" curve) does in fact dominate over 
the SM modes in the tan/3 range of interest. Note also that the combined BR's to all the 
sleptons remains under (and usually well-under) 2% even though (with ~ 150 GeV) such 
decay modes are open.0 Charged Higgs boson decays into sfermions (squarks and sleptons), 
H~ — ► qq*,£u*, are in fact heavily suppressed compared to the chargino-neutralino decay 
modes (by ~ M w ±/rriH±)', and, typically, these BR's do remain below the percent level. 

In Fig. 2, the LEP2 bound on the mass of the chargino will exclude tan/3 values above a 
certain cut-off value. The exact value of this LEP2 boundQ depends slightly on the mass of the 



electron-flavor sneutrino. For a heavy z/ e , the current LEP2 bound is m~± > 103.2 GeV [^B 



If my e < 200 GeV, this bound is weakened by only a GeV or so; however, this small change 
is enough to shift the upper limit on allowable values for tan/3 from ~23 to ~39. Also, as 
noted earlier, low values of tan/3 are excluded by LEP2 searches for a (neutral) Higgs boson. 
For this particular set of MSSM input parameters, we derive bounds of roughly tan/3 > 
2.8,2.4,2.2(3.5,2.9,2.8) for m H ± = 200, 300, 400 GeV based on the current (potential) LEP 



SM Higgs boson mass bounds given in the first (second) paper of [14]. Therefore, 3 < tan/3 < 



10, the region where charged Higgs boson signatures from SM decays are virtually absent, 
and thus the region of primary interest to us in studying chargino-neutralino decay modes, 
is not excluded. 

In this work, to avoid the enormous QCD background, only leptonic decays of the SUSY 
particles (sparticles) involved are considered. [| Two specific signal types are analyzed: events 
containing either one or three hard leptons^ accompanied by missing transverse momentum, 
p™ lss , and a reconstructed top (meaning a t- or t-quark). The top resonance is identified 
through the invariant mass of the three (at least at the parton level) jets resulting from its 
hadronic decay. (Here, we consider the rate for mis-identifying tops as very low and disregard 
any backgrounds that could arise from such mis-identification.) Tops decaying leptonically 
into a fo-jet along with a charged lepton and a neutrino are not deemed to be part of the 
signal processes, but do play a role in potentially serious backgrounds and are included in 
all simulations (as are hadronically-decaying charginos and neutralinos) . The leptons result 

5 Exceptions to this general rule are found if the stau masses are lowered to the edge of the LEP2 excluded 
region (mfi ~ 90 GeV) and m H ± <200 GeV — see Borzumati and Djouadi in j3(J • 

6 This does drop considerably if the chargino becomes near degenerate with the LSP; however, nowhere 
in the regions of parameter space we will investigate does this occur. 

7 Similar decays of neutral MSSM Higgs bosons were studied in p[ [37|. 

8 Hcrcaftcr, 'leptons' will refer to electrons and muons in general and irrespective of sign; taus and neu- 
trinos are not included. 
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from the following three-body chargino and neutralino decays: 



xt-Xi^S and (4) 

Note that the charginos and neutralinos decay directly to the LSP. Appropriate BR's for 
these decays are incorporated if necessary, but these are usually the only available decay 
modes. If the only virtual intermediate particles involved in these decays are the W ± and 
the Z°, then the leptonic BR's are the well-known leptonic BR's of the intermediate vector 
bosons (0.212 and 0.067, respectively). However, if sleptons are relatively light, they can 
also mediate these decays and significantly enhance the leptonic BR's |38| (especially those 
of the neutralinos), and thus also the rates for our prospective signals. Such light sleptons 
(with m~ oft ~ 150 GeV) are not excluded experimentally and would not be out of place in 
the light MSSM sparticle spectrum under consideration. 

In choosing the amount of missing transverse momentum required by the cuts, some 
knowledge of the chargino and neutralino mass spectrum is presumed to be available from 
independent measurements However, the charged Higgs boson mass is treated as a 



completely unknown parameter. Furthermore, due to the multiple particles leaving the 
detector unobserved, reconstruction of the Higgs boson and sparticle masses involved in the 
signal decays is not possible. Rather we here content ourselves with looking for excesses in 
the specified modes above the SM expectations. 

The analysis presented here is confined to the parton level onlyQ: jets are identified 
with the partons from which they originate and jet selection criteria are applied directly 
to the partons. Typical detector resolutions (and range limitations) are included: the 
transverse momenta of all visible particles in the final state have been smeared accord- 
ing to a Gaussian distribution, with (o"G»t)/pt) 2 = (0.6/ v /pr) 2 + (0.04) 2 for all jets and 
{(j{Pt)/pt) 2 — (0-12/y/pr) 2 + (0.01) 2 for the leptons. The missing transverse momentum has 
been evaluated from the vector sum of the jet and lepton transverse momenta after resolution 
smearing. For reference, the CTEQ4L [[13] structure function set is used, with the factoriza- 



tion scale set to Q = m t + m H ± for the Higgs boson processes and Q = m t for all others. 
Aside from using running quark masses and loop-corrected Higgs boson masses, other higher 
order corrections to the tree-level top- if 1,1 production JH] and hadronic H + — > tb decay [{45 



(which competes with our preferred SUSY decay modes) are not taken into account. The 
literature indicates that these corrections will not change the results much, though with 
small signals they should be kept in mind. 

Following values given in fl46"| , we assume a single 6-tagging efficiency of e& = 0.5 and a 
mis-tagging rate of e m i S = 0.02 (though we note the latter value may be too low since the 
study in |46|] did not include c-quarks). However, the 6-quark manifest in gg — > btH~ is 



often expected to be soft and/or near the beam pipe. Thus, = 0.5 is probably a serious 
over-estimation in this case. It is inappropriate to graft a serious 6-tagging study onto this 
parton-level analysis. A more thorough treatment will presented in the upcoming event 



9 Although the signal processes have already been incorporated into the event generators HERWIG |4Q| 
and IS A JET , the incorporation of several important background processes (see below) is still in progress 
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generator analysis |[42|| . As a simple approximation we adopt an on-off switch: if a 6-jet 
(recall this is equivalent to a bottom quark) has a pt above a certain specified value and an 
\rj\ below another specified value, we assign a 6-tagging efficiency of e b = 0.5 to it; otherwise, 
we set €b = for that soft and/or too close to the beam pipe 6-jet. Those 6-jets stemming 
from top decays are expected to almost always pass this test, so we assume e& = 0.5 for all 
such fe-jets. Now, fortuitously, it so happens that, for the particular case of one versus two 
6-jets and e b = 0.5, it does not matter how often there is one fr-jet versus how often there 
are two 6-jets fulfilling such criteria, since e b = 2e b (l — e b ) = 0.5, and so the overall 6-tagging 
efficiency for the event will remain |. This is the case for the signals we are searching for as 
well as for all the backgrounds we will discuss. 

We require one 6-tagged jet in each event. This in fact reduces both the event rates and 
the signal to background ratios for the backgrounds we will consider explicitly. However, it 
also aids in triggering and the suppression of incidental QCD backgrounds which we do not 
attempt to calculate. In addition, we do veto events with more than one 6-tagged jet. This 
does help in background reduction for all the backgrounds we consider as well as eliminating 
other bottom-rich event-types (such as gb — > ttb, gb — > tib, gg — > ttbb, etc.). We do not 
identify individual 6-quarks as tagged or untagged, rather we multiply the event by a factor 
consistent with the values given above. Again, a more technical treatment is inconsistent 
with this parton-level analysis and will come with the event generator studies. 

The one-lepton signature: + p™ lss + t 

If the charged Higgs boson is just above the top threshold, m#± ~ 200 GeV, then it is quite 
likely that the only chargino-neutralino decay channel open will be H~ — > XiXu as i s the 
case in Fig. 2. Schematically, the one lepton plus top signal would result from the reaction 
chain 

gg btH~, t bqq , H -> xIXi, Xi -> Xi^'^e , (5) 

[I = e, fi and q — d,u, s, c). The hard lepton is derived from the decay of the chargino. The 
leptonic BR of the chargino is generally not as strongly affected by a light slepton as are 
those of the non-LSP neutralinos, though modest enhancement over the expectations from 
W /± -mediated decays are possible. 

Fig. 3 gives contour plots for BR(H ± — * XiXi) with m H ± = 200 GeV, tan/5 = 4, and 
varying M 2 and fi — as noted earlier, these are the main other MSSM parameters to which 
the chargino and neutralino properties are sensitive. The region of parameter space excluded 
by the LEP2 bound on the chargino mass is indicated by the dotted curves which are, from 
bottom to top, contours for m~± = 100, 105, and 110 GeV. Note that the sensitivity to the 
exact bound here is much less than that of the tan/3 variable at the upper end of its range 
(see discussion of Fig. 2). BR's for the desired charged Higgs boson decay channel in excess of 
60% are possible in unexcluded MSSM parameter space even with this relatively low charged 
Higgs boson mass. Guided by the study of this BR, we have selected the following point in 
the MSSM parameter space for detailed simulations of the phenomenology of the one-lepton 
signature: 

M 2 = 115 GeV, n = -200 GeV, tan/? = 4 , m H ± = 200 GeV . (6) 
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At this point, relevant masses and BR's are: 



m~± = 112.61, 231.73 GeV, m~ = 59.86, 111.16, 219.43, 221.63 GeV, 
BR(H- -> xTXi) = 0.56 , BR(H~ -►&£) = 0.36 , 

BR(xr x?^) = °- 28 ( for m f ft = 150 GeV ) • ( 7 ) 

SM backgrounds to such events come from top pair production and single top production: 

gg, qq->ti , t -> bqq' , i -> 6£"z/, , (8) 
5-5-, gg -> tW" , t -> &g<f , W~ -> f z/, , (9) 

where the initial pair in (|j) does not come from an on-shell top decay. At the LHC, 
approximately 0.1 billion ti events will be produced for every 100 fb -1 of integrated luminos- 
ity; whereas, the corresponding number of top- if ± events is only several thousand — note 
that (|) is suppressed relative to (D) by ~ a em , meaning that rates for both backgrounds are 
larger than that of the would-be signal. Finally, the bottom-top decay of the charged Higgs 
boson may yet have an appreciable BR even when chargino-neutralino decay modes are very 
important. Such C H~ — ► bt ' events; i.e., 

gg — > btH~ , t — > bqq , H~~ — > bt , i — > bl~v t 
or gg -> btH~ , t -> 6£ + z/, , F~ -»• 6f , f -> , (10) 

might also pass our signal cuts, though these are not designed to optimize the selection of 
H~ — > bt events which, for instance, have four 6-jets manifest in the decay chains of (|l0|). 
whereas only one tagged 6-jet is permitted by our selection criteria. To cut out additional 
QCD backgrounds, such as the radiation of hard gluons from the afore- mentioned SM back- 
grounds (or MSSM gluino pair production followed by cascade decays), we will put a 4-jet 
cap[^] on the number of jets we allow in any event. Adding the four 6-jets above and two 
distinct (neglecting the rare case of jet mis-identification) untagged jets we will require to 
reconstruct a hadronically-decaying yields six jets, meaning most H~ — > bt events will 
also be lost to the 4-jet cut.Q Combined with the compulsory single tagged 6-jet, this also 
implies that the surviving H — > bt events will have at least one and at most two 6-jets 
passing our on-off switch criteria, and so the 6-tagging efficiency factor for these surviving 
events will again be |. 

In our simulation, we have adopted the following acceptance and selection cuts: 

10 Of course, the jet number in an event can be affected, for example, by the merging of showers from 
different initial partons or hadrons lying too close to the beam pipe. A full event-generator analysis should 
yield a more accurate estimate of the fraction of the time an expected jet is not seen than the present 
parton-level analysis; so herein we will neglect such additional backgrounds. 

11 Because of this, associated production, gg, qq — > W + H~ , with the W + providing the hard lepton and the 
top coming from H~ — > bt, might be a mimic of comparable size to (||), our designated reaction chain, (even 
though the W^H^ production cross-section at this MSSM point is down by roughly an order of magnitude 
from top-H^ production) if a higher fraction of such events survive the cuts. However, since we will show 
explicitly that the reaction chains in ( ]To| ) have a negligible effect, it is clear that this alternative reaction (or 
the even more suppressed charged Higgs pair production reaction chains) will also be unimportant. 
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1. Jets and leptons are retained if they satisfy the following requirements: pj? > 25 GeV, 
\r/£j\ < 2 and AR^j/jj = \J^Vej/jj + ~^$fj/jj > where j represents both 6-tagged 
and untagged jets. To pass this first cut, an event must have one and only one lepton 
that satisfies the above criteria and no more than four jets (irrespective of whether or 
not the jets are 6-tagged) fulfilling the requirements. 

2. We require that pf ss > 80 GeV. 

3. We demand that two untagged jets reproduce an invariant mass around M w ±: 

\M q qi — M w ± | < 10 GeV. Therefore, to take into account the possibility of mis-tagging 
a jet as a 6-jet, we will multiply the 6-tagging factor to be applied at the end of the 
series of cuts by (1 — 2e mis + e^ is ) = 0.96 . 

4. We combine these two light-quark-jets with a 6-jet and demand that at least one such 
resulting 3-jet invariant mass be in the vicinity of m t : {M^gi — m t \ < 25 GeV. 

5. Recognizing the difference in the number and type of particles leaving the detector 
undetected in the signal events (which have two LSP's as well as neutrinos) and back- 
ground events (which have only neutrinos - no LSP's), we construct a variable to 

exploit this distinction, demanding that j p m iss+p £j > 0.2. 

6. Finally, we apply a veto on a leptonically decaying t- or t-quark. If more than three jets 
(that is, jets in addition to the three jets already assigned to a hadronically-decaying 
top in 3. and 4.) are present, then an invariant mass denoted by M&&^ is formed from 
each extraneous jet's four-momentum and those of the hard lepton and the missing 
momentum. The missing momentum is assumed to be solely due to a massless neutrino 
whose longitudinal momentum is reconstructed following the technique outlined in the 



first paper of |31[] . This assignment is quite reasonable for SM double- and single-top 
events,^ so that the former can be eliminated with a cut of \M^ ul —m t \ > 25 GeV. The 
assignment is of course not at all reasonable for decays of the charged Higgs boson, and 
simulations confirm that a far smaller fraction of these events are lost in comparison 
to the percentage of single- and double-top events weeded out. 

Results for this series of cuts are given in Tab. (EJ With an integrated luminosity of 
100 fb -1 , about 265 signal events per year of run time survive. However, approximately 
50,000 background events also survive. Even if one considers more favorable points in the 
MSSM parameter space, it is very difficult to enhance the signal cross-section significantly. 
Thus, if this one-lepton channel is to be useful in searching for charged Higgs boson at the 
LHC, far better cuts than those designed here will need to be devised. 



12 Since neutrinos that may produced in decays of B- mesons or further on down the decay chains inside 
the tagged and untagged &-jets are generally fairly soft. 
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oOo. 


OD. 




i on 


< 4 jets with > 25 GeV, 










\r)j\ < 2, AR J:j > 0.4 


153. 


16. 


.033 


.042 


pf ss > 80 GeV 


28. 


3.53 


.0049 


.025 


|M 5r - M w ± \ < 10 GeV 


27. 


3.38 


.0040 


.019 


\M bq q>-m t < 25 GeV 


25. 


2.93 


.0037 


.018 


1 miss £ 1 

|p™ iss +p*,| ^ u ' z 


18. 


2.65 


.0029 


.017 


|M 6 ^ -m t | > 25 GeV 


1.99 


1.49 


.0005 


.014 



Table 1: Production and decay rates (in picobarns) for one-lepton signal and backgrounds, after 
the implementation of successive cuts. Rates already include the Higgs boson decay BR's, whereas 
the common 6-taggingx mis-tagging factor of ^ x 0.96 is not included. Also omitted are the BR 
for one top to decay hadronically (0.699) and the leptonic BR's, 0.212 for the top decays and 
BR(x7 ~~ > Xi^~Ve) for the signal, as well as the consequent combinatorial factor of 2. 



The three-lepton signature: t^£ £ + + p™ lss + t 

For larger charged Higgs boson masses, other chargino-neutralino decay modes besides 
H~ — > XiXi may well be open and have sizable BR's (as shown in Fig. 2). The heavier 
neutralinos may then decay leptonically (f|) and, together with the chargino, produce three 
hard leptons (as first discussed in P4"| ). Recall that in our notation T stands for either 



electrons or muons. For the signal, two leptons with opposite signs must be of the same 
flavor; the third lepton may also be of the same flavor or of the other flavor. The expected 
reaction chain for the signal is 

gg^btH', t^bqq', H~ - xiX°,3 , Xx - Xx^ , xh^Xx^^. (11) 

Fig. 4 gives contour plots for BR(if ± — > xfxV) ( on top) and BR(if ± — > xtx.3) ( on bottom), 
with mii± = 300 GeV, tan/3 = 4, and again varying M 2 and [i. Limits of the regions excluded 
by LEP2 are again marked by dotted lines. BR's above 20% are found at viable points in the 
parameter space. (The diagonal discontinuities seen in the upper right corners of the fi < 
plots and the upper left corners of the fi > plots are due to a 'level-crossing' where the 
masses of x° an d X3 become degenerate — thus the identities of these two neutralinos are 
effectively interchanged as one of these diagonal lines is crossed.) From the study of these 
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BR's the following point in the MSSM parameter space was chosen for the simulation study: 



M 2 = 200 GeV, p = -120 GeV, tan/3 = 4, m H ± = 300 GeV, mf ft = 150 GeV . (12) 
At this point, relevant masses and BR's are: 

m~± = 116.85, 231.48 GeV, m~„ = 87.93, 122.01, 140.29, 230.76 GeV, 

<M,2 Al_4 

BR(H~ - xrX2( 3) ) = 0.18(0.03) , BR(H~ - 6f) = 0.63 , 

BR(X2( 3 ) - Xi^O = 0.33(0.02) , BR(xr - Xi^P,) = 0.24 . (13) 
Additional key variables to be aware of are m~- — m~o and m~o — m~o. These are not 

Xi Xi X 2 ,3 Xi 

so large here and this softens both the lepton spectra (the leptons coming from chargino 
and neutralino decays have on average lower transverse momenta than those coming from 
gauge boson decays) and that of p™ lss . Another point worthy of mention is that the MSSM 
parameter point ([12]), gives rrih = 105.5 GeV if soft stop masses are set to 1 TeV and A t = 0. 
With a projected LEP2 reach of E cm ~ 208 GeV, this would yield a Higgsstrahlung cross- 
section that should be observable. However, if, for instance, A t is raised to 2 TeV, then 
rrih = 118.9 GeV and on-shell Higgsstrahlung is kinematically forbidden. While soft SUSY- 
breaking parameters such as A t can have a strong impact on rrih (as well as a possibly 
significant impact on the b — > 57 rates, as will be discussed later), they have very little effect 
on rriH± (as noted earlier). Thus some care must be taken that all relevant free parameters 
in the model are adequately explored so as to not neglect allowable MSSM parameter sets. 

The dominant SM backgrounds are again those involving double- and single-top produc- 
tion and decay, this time accompanied by an additional lepton-antilepton pair (electrons or 
muons) produced in the 'off-shell decay' of a neutral gauge boson {V = j/Z)^ 

gg, qq -»■ tiV* , t^bqq' , t -> U~v t , V* -> t ~E + , (14) 
and gg, qq -> tbW~V* , t -> bqq' , W~ -> rp, , V* -> . (15) 

The set of cuts applied is similar that employed for the one-lepton signal analysis: 

1. Jets and leptons are retained if they satisfy the following requirements: p l T > 10 GeV, 
Pr > 25 GeV, \i]ej\ < 2 and ARij/jj > 0.4. The signal rate is sensitive to the pt 
threshold for the leptons — lowering the threshold will enhance the signal survival 
rate more than that of the backgrounds. The value chosen here is reflective of the 
capabilities of the ATLAS detector fl6]| . As with the cuts for the one-lepton signal, for 
an event to pass this first cut it is compelled to have no more than four jets fulfilling 
the requirements. And in this case we of course demand that exactly three leptons also 
satisfy the criteria. 

2. We require that pf ss > 25 GeV. 



13 To assess the ttV* background, the code originally developed in Q was adapted to allow for an off-shell 
;;auge boson. 
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3. As before, we impose \M q qi — M w ± \ < 10 GeV. 

4. We also again impose \Mbqqi — m t \ < 25 GeV. 

5. Given the rather low missing momenta involved in both signal and backgrounds, we 
find that the variable | p m iss+p f used in the one-lepton analysis is no longer a suitable 
discriminant and do not include it in the cuts. 

6. Lastly, we apply a Z°-veto, \M^-^+ — Mz\ > 10 GeV. This is to eliminate the SM 
backgrounds where the gauge boson is on- or nearly on-shell. Note that, for the signal, 
m~o — m~o <C Mz, and thus few signal events are lost here. 

Results for this series of cuts are given in Tab. |2|. Unlike in the case of the one-lepton 
signature, after cuts the three-lepton signal rate can be made competitive with the back- 
ground rates. However, the total number of signal events is low. If one multiplies the 
final row of numbers of Tab. Q by the 6-taggingx mis-tagging factor of \ x 0.96 and by the 
leptonic W and Z° BR's (for first two columns) or the leptonic branching ratios from 
Eqs. (|1^) (for last two columns), one finds for one year's running (100fb _1 ) the ratio, 
signal events : background events = 7 : 20. Note that the enhanced leptonic BR of x° 
due to the light slepton intermediate state is very significant. 

A more stringent cut can be applied on the invariant mass of the opposite-sign lepton 
pair if we take advantage of the fact that M^-£+ < m x ( > 3 ~ m xl < Mz- This entails the 
possession of some information about the masses of the lowest-lying sparticle states. As 
noted earlier, it is quite likely that such information will be available to those analyzing 
the real experimental data in search of a signal for the charged Higgs boson. (Note that it 
should also be possible to tune the M1-1+ limit to help optimize any observed signal even 
with incomplete information about the sparticle masses.) Seeing that the X1X3 contribution 
is not very important, we could impose the cut M^- t + < m~o — m~o, with the inclusion of 
which one finds for one year's running the ratio, signal events : background events = 5:5. 
Note that with the parameter set (|T2"D and the consequent sparticle mass spectrum (0), cut 
6. of Tab.||] will be completely subsumed by this new cut. While this should be the case in 
general, it is safer to separately apply cut 6. and this new cut since the choice of the former 
is not parameter-space dependent and the cut-off for the later might rise above Mz — 10 GeV 
in exceptional cases. 

Another additional cut which has some dependence on the mass spectrum of the charginos 
and neutralinos can be applied by defining Mt(3£) = \j2p^p^ lss (l — cosA</>), where p^ is 
the transverse momentum of the three-lepton system and A</> is the azimuthal separation 
between p^ and p™ lss . Fig. 5 illustrates how well this variable distinguishes our signal from 
the backgrounds. For the former, the M T (3£) distribution dies at m#± — 2m~o, which is 
~ 123 GeV for the point (p~2"D; whereas, for the latter, it can stretch far beyond this value. 
Demanding Mt(3£) < 100 GeV yields the ratio, signal events : background events =7:9 
(5 : 2), if applied on top of the cuts in Tab. |2] (or in conjunction with strengthening the 
\M t - l+ - M z \ > 10 GeV cut in Tab. g to M e - i+ < m~o - m~o). 
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ttv* 


tbW'V* 


H- -> Xixl 


H - Xixl 


No cuts 


698. 


111. 


43. 


7.0 


3 lept 


ons each with p l T > 10 GeV, 
Vt\<2, AR(£,j) > 0.4 


317. 


49.7 


8.7 


2.2 


< 4 jets with p{ > 25 GeV, 
\Vj\<2, AR(j,j)>0A 


161. 


30.3 


1.75 


.43 


pf ss > 25 GeV 


133. 


21.1 


1.67 


.42 




M q g> - M w ± 


< 10 GeV 


126. 


20.9 


1.46 


.39 




\M bqq i - m t 


< 25 GeV 


110. 


11.3 


1.41 


.33 




M e - e+ - M z 


> 10 GeV 


17. 


5.03 


1.38 


.32 



Table 2: Production and decay rates (in femtobarns) for three-lepton signal and backgrounds, after 
the implementation of successive cuts. Rates already include the Higgs boson decay BR's, whereas 
a common hadronic BR for the decaying top (0.699) and a common 6-taggingx mis-tagging factor 
of ^x0.96 have been omitted. The SM backgrounds should also be multiplied by 2 x 0.212 x 0.066 
(accounting for the and Z° leptonic BR's), while the signal rates should be multiplied by 
2 x BR(xr -» T&tTDt) x BR(x° (3) Xi^ + )- The tiV* and tbW'V* cross-sections (V = j,Z) 
are expressible in terms of the Z° decay rates since, at the end of the series of cuts, the Z° — » £~£ + 
contribution is numerically dominant over the one from 7* — ► £~£ + . Also note that, since rri£ is set 
to zero, a hard cut of M e -g+ > 10 GeV is necessary to avoid the 7* — ► £~£ + singularity — this is 
included in the "No cuts" rates for tbW~V* and ttV*. 



Bolstered somewhat by this result, it is reasonable to consider still higher charged Higgs 
boson masses. Now competing factors come into play. One the one hand, the three- 
lepton chargino-neutralino decay channels remain large — with m H ± = 400 GeV and all 
other MSSM parameters the same as in point (|1~2|), one has BR(H~ — > X1X2) = 0.14 
and BR(if ~ — > X1X3) = 0.03 — and, on the positive side, the lepton spectrum becomes 
harder. On the other hand, the top-H ± production rate drops precipitously, by more than 
a factor of 2 when m#± is increased from 300 GeV to 400 GeV. Unfortunately the lat- 
ter negative effect dominates: although a greater percentage of the signal events survive 
the cuts, one starts with a production cross-section that is just too small. Applying the 
same cuts as in Tab.|| for m#± = 400 GeV yields 0.55 fb and 0.20 fb for the H~ — ► X1X2 
and H — > X1X3 channels, respectively. Adopting the same 6-tagging efficiencies and lep- 
tonic BR's as for the m H ± = 300 GeV case now yields (for m H ± = 400 GeV) the ratio, 
signal events : background events = 3 : 20. In addition, the larger m#± means the extra 
M T (3£) cut must be weakened — requiring M T (3£) < 120 GeV along with \Mn- t + — M z \ > 
10 GeV leads to the ratio, signal events : background events = 3 : 11. Also, since the 
H — > xi~X3 decay modes acconts from more of the signal now, if the Mi-£+ < m~o — m~o 
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cut is implemented, more signal will be lost, while a weaker M^-i+ < m~o —m~o cut is much 

A3 K% 

less effective at cutting away backgrounds. Thus, both the total three-lepton signal event 
rate and its statistical significance decline as m#± is raised from 300 GeV to 400 GeV. To 
this though must be added the caveat that, at m H ± = 400 GeV, decay modes including ei- 
ther the heaviest chargino or the heaviest neutralino are significant (as seen from Fig. 2). We 
have neglected these, and so our results may be viewed as conservative. But it is nonetheless 
evident that m#± = 400 GeV is near the kinematical limit beyond which there is too lit- 
tle top-H^ production cross-section at the LHC to exploit through the chargino-neutralino 
decay channels. 

In contrast, for m#± <400 GeV, tan/? values up to ~ 10 can be scanned (for a significant 
portion of the possible values of the other MSSM input parameters) using the trilepton plus 
top signature from H ± — > xtx.23 decays. In fact, over the range 3 < tan/3 <10, as tan/3 
gets larger, the enhancement (or suppression) of the Higgs decay rates is compensated by an 
opposite effect in the production rate. This is understandable since the BR's for H ± — > xtx^ 
strengthen as the H — > bi decay width weakens (as can be seen from an examination of 
Fig. 2), and the latter is proportional to the same coupling (|TJ) as the top- if ± production 
modes. Beyond tan/3 ~ 10, the BR(if ± — > xf x^'s start falling below the tan/3 = 4 values 
simulated in the preceding numerical analyses. In fact, due to the strengthening of other 
alternative decay modes (such as H — > t~P t ), these BR's fall even quicker than the top- H ± 
production rate increases. 

Since the event rate imposes a discovery limit (< 400 GeV) on the mass of the charged 
Higgs bosons, it is important to determine if all or part of the region of parameter space 
where the signal seems observable is excluded by constraints coming from b — > 37 decays. 
At present however, we are unaware of any clean answer to this question. To circumvent 
the b —>■ 57 constraint with such modest charged Higgs boson masses, there must be suf- 
ficient cancellation between the contributions from the top-charged Higgs boson loop and 
the chargino-stop loop (we assume contributions from gluino-squark loops [^] are negligible, 
though this need not always be true). Seeing that our signal stems from charged Higgs boson 
decays including a chargino (this prefers, among other things, lower values for |/z|), it follows 
that the main freedom left to adjust the b — > S7 rates lies in the stop sector.^ The MSSM 
stop input parameters must be large enough to boost the light Higgs boson mass above the 
LEP2 bounds (for the sample point we have used for illustration, this would translate into 
stop masses in the ~400-600 GeV ballpark), yet low enough to ensure sufficient cancellation 
among the MSSM b — > S7 diagrams (see |TD|, Fig. 19, for a rough idea). Studies addressing 
this have been done in the more restrictive mSUGRA scenario (see for example Goto and 



Okada of ||13|1 ; note though that this is a leading-order calculation, and that more recent 
next-to-leading order studies suggest there may be significant modifications to these results 
- see Ciuchini et al. and Bobeth et al. of JT3I), and point toward a severe curtailment of the 



range of open parameter space for our signal. However, in mSUGRA the stop parameters 



14 Here we are operating under the assumption that inter-generational squark mixing is negligible, anal- 
ogous to the case for the third generation of the CKM matrix. If squark mixings are made arbitrary, the 
b — > sj constraint can almost certainly be evaded; however, constraints from other flavor-changing neutral 
current processes need to be considered as well. 



13 



are related to those of the Higgs bosons and the Higgsinos. If a more general MSSM scenario 
is assumed these relationships will disappear and the restriction from b — > may well be 
significantly relaxed. 

One final issue to consider is the distinguishability of the charged Higgs bosons from the 
neutral H and A. Throughout much of the parameter space, these MSSM Higgs bosons all 
have very similar masses. Background to the charged Higgs boson signal could come from 
gg,qq -> ttH,tiA |g g[J or qq' -> tbH,tbA production 0, where the neutral MSSM Higgs 
boson then decays into chargino or neutralino pairs. Fortunately, for the moderate tan/? 
values we are interested in here, the tiH and tiA production rates at the LHC are about an 
order of magnitude lower than the top- if ± production rate. (See also paper #2 of fl32f .) 

In summary, chargino-neutralino decays of heavy (i.e., with raa± Zm t ) MSSM charged 
Higgs bosons have some potential to aid in the detection of such H ± scalars at the LHC, 
especially in the intermediate tan/3 region, 3 < tan/5 < 10, which is inaccessible via SM de- 
cays.0 Among the possible chargino-neutralino combinations, the most promising are the 
H~ -> XiXliXiXl decays followed in turn by the decays Xi -> &~v t and x°, 3 
yielding a three-lepton final state. (The decay H~ — ► XiXa leading to a one-lepton final 
state was also studied, but found to be overwhelmed by the SM backgrounds.) Since the 
charged Higgs boson is dominantly produced in association with a top quark, the main sig- 
nature to look for is three hard, isolated leptons plus significant p™ ss and a reconstructed 
top quark. Unfortunately, after the cuts utilized in this simulation, the surviving signal 
event rates are small — only around a handful of events per year. Furthermore, these rates 
are sensitive to several MSSM parameters, including M 2 , /x, and (via the leptonic BR's of 
X.2 3) m |° ft - However, given the paucity of handles which can be used to study charged Higgs 
boson production at the LHC for intermediate tan/3 values, further investigation of even such 
weak signals in the more realistic environment of a full event generator, including parton 
shower effects and hadronization, would be clarifying and beneficial. Such a study is well 



underway (using the HERWIG [40] and ISAJET |4T[ event generators) and we plan to report 



the results from these studies in the near future |4"2" | 
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Figure 1: BR's of the MSSM charged Higgs boson into chargino-neutralino pairs (summing all 
such channels), in the (rattan/?) plane, with M 2 = 200 GeV and \x = — 120 GeV. One-loop 
formulas as found in [33], |41| are used to relate m#± to tua- Other MSSM input parameters 
are: m| oft = 1 TeV, A t = 2 TeV and m| oft = 300 GeV. 
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Figure 2: BR's of the MSSM charged Higgs boson as a function of tan/3 for m H ± = 200, 300, 
and 400 GeV, again with M 2 = 200 GeV and /i = -120 GeV. Other MSSM input parameters 
are also as in Fig. 1, except that here m| oft = 150 GeV. 




Figure 3: BR(H + -> XiXi) in the (/i, M 2 ) plane for tan/3 = 4 and m H ± = 200 GeV. Other 
MSSM input parameters are as in Fig. 1. The dotted lines for m ± indicate the reach of 

Ai 

LEP2. As expected, setting m~ oft = 150 GeV or m~ oft = 300 GeV does not affect the results. 




Figure 4: BR(H + — > Xi^x") ( on to P) an d BR(i/+ — ► X1X3) ( on bottom) in the (fx, M 2 ) plane 
for tan/5 = 4 and m#± = 300 GeV. Other MSSM input parameters are as in Fig. 1, except 
that m| oft = 150 GeV. The dotted lines for m + indicate the reach of LEP2. 



BR(H + -y 


h 


l yoh\ 

2 \ K 
0.15 X 


"Olo2 








0.1 0.2 


Hj 


_ / / 


/ 


0.05' 


\ \ w 


^25 


/ ( 


).02 








I0.01 










m z ,(in GeV) :> 


\I10 
1B5 
10G7 








---i--t~T" \ 


r'V 




'i IN I > 





250 
230 
210 

190 
1 70 
1 50 
1 30 
1 10 l 



ho3i 




\ \\ 


30.2 


\ 












\ 0.02\0.01" 
005 \ \ 


^7 0.15/ 






I 1 1 1 1/ 


0.1/ 

1ii 


1 1 1 1 1 1I1 1 1 1 1 1 



-200 -175 -150 -125 -100 125 150 175 200 225 250 



tanfi = 4 



[l ( GeV) m H+ = 300 GeV 



250 




D 2i \K 005 \ 
-K\ 0.02 


\0.01 


— \\o.t\ 






\A m 7 +(in GeV) > 110 
■••\\ Vl05 






I i i i i i K 


I \ i\l \\ i I i i i 



200 -175 -150 -125 -100 125 150 175 200 225 250 



\x (GeV) 



Figure 5: Normalized differential distributions of the three-lepton system transverse mass, 
M T (3£) (as defined in the text) for: H + — > X1X2 (solid: m H ± = 300 GeV; dashed: m H ± = 
400 GeV), H + -> xtxl (fine-dotted: m H ± = 300 GeV; dot-dashed: m H ± = 400 GeV), ttV* 
(long-dashed), and tbW-V* (dotted). Here M 2 = 200 GeV, n = -120 GeV and tan/3 = 4. 
Other MSSM input parameters are as in Fig. 1. Normalizations are as after the cuts in 
Tab.|2]; i.e., the leptonic BR's and the fe-taggingx mis-tagging efficiency are not included. 




